Analysis of somatic mutations in whole blood from 200,618 individuals identifies pervasive positive selection and novel drivers of clonal hematopoiesis.

Human aging is marked by the emergence of a tapestry of clonal expansions in dividing tissues, particularly evident in blood as clonal hematopoiesis (CH). CH, linked to cancer risk and aging-related phenotypes, often stems from somatic mutations in a set of established genes. However, the majority of clones lack known drivers. Here we infer gene-level positive selection in whole blood exomes from 200,618 individuals in UK Biobank. We identify 17 additional genes, ZBTB33, ZNF318, ZNF234, SPRED2, SH2B3, SRCAP, SIK3, SRSF1, CHEK2, CCDC115, CCL22, BAX, YLPM1, MYD88, MTA2, MAGEC3 and IGLL5, under positive selection at a population level, and validate this selection pattern in 10,837 whole genomes from single-cell-derived hematopoietic colonies. Clones with mutations in these genes grow in frequency and size with age, comparable to classical CH drivers. They correlate with heightened risk of infection, death and hematological malignancy, highlighting the significance of these additional genes in the aging process.

Reconstructing phylogenetic trees from genome-wide somatic mutations in clonal samples.

Phylogenetic trees are a powerful means to display the evolutionary history of species, pathogens and, more recently, individual cells of the human body. Whole-genome sequencing of laser capture microdissections or expanded stem cells has allowed the discovery of somatic mutations in clones, which can be used as natural barcodes to reconstruct the developmental history of individual cells. Here we describe Sequoia, our pipeline to reconstruct lineage trees from clones of normal cells. Candidate somatic mutations are called against the human reference genome and filtered to exclude germline mutations and artifactual variants. These filtered somatic mutations form the basis for phylogeny reconstruction using a maximum parsimony framework. Lastly, we use a maximum likelihood framework to explicitly map mutations to branches in the phylogenetic tree. The resulting phylogenies can then serve as a basis for many subsequent analyses, including investigating embryonic development, tissue dynamics in health and disease, and mutational signatures. Sequoia can be readily applied to any clonal somatic mutation dataset, including single-cell DNA sequencing datasets, using the commands and scripts provided. Moreover, Sequoia is highly flexible and can be easily customized. Typically, the runtime of the core script ranges from minutes to an hour for datasets with a moderate number (50,000-150,000) of variants. Competent bioinformatic skills, including in-depth knowledge of the R programming language, are required. A high-performance computing cluster (one that is capable of running mutation-calling algorithms and other aspects of the analysis at scale) is also required, especially if handling large datasets.

Inherited blood cancer predisposition through altered transcription elongation.

Despite advances in defining diverse somatic mutations that cause myeloid malignancies, a significant heritable component for these cancers remains largely unexplained. Here, we perform rare variant association studies in a large population cohort to identify inherited predisposition genes for these blood cancers. CTR9, which encodes a key component of the PAF1 transcription elongation complex, is among the significant genes identified. The risk variants found in the cases cause loss of function and result in a ∼10-fold increased odds of acquiring a myeloid malignancy. Partial CTR9 loss of function expands human hematopoietic stem cells (HSCs) by increased super elongation complex-mediated transcriptional activity, which thereby increases the expression of key regulators of HSC self-renewal. By following up on insights from a human genetic study examining inherited predisposition to the myeloid malignancies, we define a previously unknown antagonistic interaction between the PAF1 and super elongation complexes. These insights could enable targeted approaches for blood cancer prevention.

Inherited polygenic effects on common hematological traits influence clonal selection on JAK2V617F and the development of myeloproliferative neoplasms.

Myeloproliferative neoplasms (MPNs) are chronic cancers characterized by overproduction of mature blood cells. Their causative somatic mutations, for example, JAK2V617F, are common in the population, yet only a minority of carriers develop MPN. Here we show that the inherited polygenic loci that underlie common hematological traits influence JAK2V617F clonal expansion. We identify polygenic risk scores (PGSs) for monocyte count and plateletcrit as new risk factors for JAK2V617F positivity. PGSs for several hematological traits influenced the risk of different MPN subtypes, with low PGSs for two platelet traits also showing protective effects in JAK2V617F carriers, making them two to three times less likely to have essential thrombocythemia than carriers with high PGSs. We observed that extreme hematological PGSs may contribute to an MPN diagnosis in the absence of somatic driver mutations. Our study showcases how polygenic backgrounds underlying common hematological traits influence both clonal selection on somatic mutations and the subsequent phenotype of cancer.

Diagnosis and evaluation of prognosis of myelofibrosis: a British Society for Haematology guideline

This document represents an update of the British Society for Haematology (BSH) guideline on myelofibrosis (MF) first published in 2012 and updated in 2015.1 This guideline aims to provide healthcare professionals with clear guidance on the diagnosis and prognostic evaluation of primary my-elofibrosis (PMF), as well as post-polycythaemia vera myelo-fibrosis (post-PV MF) and post-essential thrombocythaemia myelofibrosis (post-ET MF). A section on prefibrotic MF is also included. A separate BSH Guideline covers the manage-ment of MF and is published alongside this guideline.

The management of myelofibrosis: a British Society for Haematology guideline

This document represents an update of the British Society for Haematology guideline on Myelofibrosis first published in 2012 and updated in 2015 These guidelines aim to pro-vide healthcare professionals with clear guidance on stratified management for primary myelofibrosis (PMF), as well as post-polycythaemia myelofibrosis (post-PV MF) and postessential thrombocythaemia myelofibrosis (post-ET MF). A separate BSH guideline covers the diagnosis and prognostic evaluation of myelofibrosis and is published alongside this guideline

Clonal selection of hematopoietic stem cells after gene therapy for sickle cell disease.

Gene therapy (GT) provides a potentially curative treatment option for patients with sickle cell disease (SCD); however, the occurrence of myeloid malignancies in GT clinical trials has prompted concern, with several postulated mechanisms. Here, we used whole-genome sequencing to track hematopoietic stem cells (HSCs) from six patients with SCD at pre- and post-GT time points to map the somatic mutation and clonal landscape of gene-modified and unmodified HSCs. Pre-GT, phylogenetic trees were highly polyclonal and mutation burdens per cell were elevated in some, but not all, patients. Post-GT, no clonal expansions were identified among gene-modified or unmodified cells; however, an increased frequency of potential driver mutations associated with myeloid neoplasms or clonal hematopoiesis (DNMT3A- and EZH2-mutated clones in particular) was observed in both genetically modified and unmodified cells, suggesting positive selection of mutant clones during GT. This work sheds light on HSC clonal dynamics and the mutational landscape after GT in SCD, highlighting the enhanced fitness of some HSCs harboring pre-existing driver mutations. Future studies should define the long-term fate of mutant clones, including any contribution to expansions associated with myeloid neoplasms.

Convergent somatic evolution commences in utero in a germline ribosomopathy.

Clonal tracking of cells using somatic mutations permits exploration of clonal dynamics in human disease. Here, we perform whole genome sequencing of 323 haematopoietic colonies from 10 individuals with the inherited ribosomopathy Shwachman-Diamond syndrome to reconstruct haematopoietic phylogenies. In ~30% of colonies, we identify mutually exclusive mutations in TP53, EIF6, RPL5, RPL22, PRPF8, plus chromosome 7 and 15 aberrations that increase SBDS and EFL1 gene dosage, respectively. Target gene mutations commence in utero, resulting in a profusion of clonal expansions, with only a few haematopoietic stem cell lineages (mean 8, range 1-24) contributing ~50% of haematopoietic colonies across 8 individuals (range 4-100% clonality) by young adulthood. Rapid clonal expansion during disease transformation is associated with biallelic TP53 mutations and increased mutation burden. Our study highlights how convergent somatic mutation of the p53-dependent nucleolar surveillance pathway offsets the deleterious effects of germline ribosomopathy but increases opportunity for TP53-mutated cancer evolution.

Ruxolitinib Versus Best Available Therapy for Polycythemia Vera Intolerant or Resistant to Hydroxycarbamide in a Randomized Trial.

PURPOSE: Polycythemia vera (PV) is characterized by JAK/STAT activation, thrombotic/hemorrhagic events, systemic symptoms, and disease transformation. In high-risk PV, ruxolitinib controls blood counts and improves symptoms. PATIENTS AND METHODS: MAJIC-PV is a randomized phase II trial of ruxolitinib versus best available therapy (BAT) in patients resistant/intolerant to hydroxycarbamide (HC-INT/RES). Primary outcome was complete response (CR) within 1 year. Secondary outcomes included duration of response, event-free survival (EFS), symptom, and molecular response. RESULTS: One hundred eighty patients were randomly assigned. CR was achieved in 40 (43%) patients on ruxolitinib versus 23 (26%) on BAT (odds ratio, 2.12; 90% CI, 1.25 to 3.60; P = .02). Duration of CR was superior for ruxolitinib (hazard ratio [HR], 0.38; 95% CI, 0.24 to 0.61; P < .001). Symptom responses were better with ruxolitinib and durable. EFS (major thrombosis, hemorrhage, transformation, and death) was superior for patients attaining CR within 1 year (HR, 0.41; 95% CI, 0.21 to 0.78; P = .01); and those on ruxolitinib (HR, 0.58; 95% CI, 0.35 to 0.94; P = .03). Serial analysis of JAK2V617F variant allele fraction revealed molecular response was more frequent with ruxolitinib and was associated with improved outcomes (progression-free survival [PFS] P = .001, EFS P = .001, overall survival P = .01) and clearance of JAK2V617F stem/progenitor cells. ASXL1 mutations predicted for adverse EFS (HR, 3.02; 95% CI, 1.47 to 6.17; P = .003). The safety profile of ruxolitinib was as previously reported. CONCLUSION: The MAJIC-PV study demonstrates ruxolitinib treatment benefits HC-INT/RES PV patients with superior CR, and EFS as well as molecular response; importantly also demonstrating for the first time, to our knowledge, that molecular response is linked to EFS, PFS, and OS.

Clonal dynamics of haematopoiesis across the human lifespan.

Age-related change in human haematopoiesis causes reduced regenerative capacity1, cytopenias2, immune dysfunction3 and increased risk of blood cancer4-6, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at constant rates throughout life, affecting more genes than identified in blood cancers. Loss of the Y chromosome conferred selective benefits in males. Simulations of haematopoiesis, with constant stem cell population size and constant acquisition of driver mutations conferring moderate fitness benefits, entirely explained the abrupt change in clonal structure in the elderly. Rapidly decreasing clonal diversity is a universal feature of haematopoiesis in aged humans, underpinned by pervasive positive selection acting on many more genes than currently identified.

The longitudinal dynamics and natural history of clonal haematopoiesis.

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis1-4. The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited5,6. Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2P95H). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection.

Bayesian networks elucidate complex genomic landscapes in cancer.

Bayesian networks (BNs) are disciplined, explainable Artificial Intelligence models that can describe structured joint probability spaces. In the context of understanding complex relations between a number of variables in biological settings, they can be constructed from observed data and can provide a guiding, graphical tool in exploring such relations. Here we propose BNs for elucidating the relations between driver events in large cancer genomic datasets. We present a methodology that is specifically tailored to biologists and clinicians as they are the main producers of such datasets. We achieve this by using an optimal BN learning algorithm based on well established likelihood functions and by utilising just two tuning parameters, both of which are easy to set and have intuitive readings. To enhance value to clinicians, we introduce (a) the use of heatmaps for families in each network, and (b) visualising pairwise co-occurrence statistics on the network. For binary data, an optional step of fitting logic gates can be employed. We show how our methodology enhances pairwise testing and how biologists and clinicians can use BNs for discussing the main relations among driver events in large genomic cohorts. We demonstrate the utility of our methodology by applying it to 5 cancer datasets revealing complex genomic landscapes. Our networks identify central patterns in all datasets including a central 4-way mutual exclusivity between HDR, t(4,14), t(11,14) and t(14,16) in myeloma, and a 3-way mutual exclusivity of three major players: CALR, JAK2 and MPL, in myeloproliferative neoplasms. These analyses demonstrate that our methodology can play a central role in the study of large genomic cancer datasets.

Life histories of myeloproliferative neoplasms inferred from phylogenies.

Mutations in cancer-associated genes drive tumour outgrowth, but our knowledge of the timing of driver mutations and subsequent clonal dynamics is limited1-3. Here, using whole-genome sequencing of 1,013 clonal haematopoietic colonies from 12 patients with myeloproliferative neoplasms, we identified 580,133 somatic mutations to reconstruct haematopoietic phylogenies and determine clonal histories. Driver mutations were estimated to occur early in life, including the in utero period. JAK2V617F was estimated to have been acquired by 33 weeks of gestation to 10.8 years of age in 5 patients in whom JAK2V617F was the first event. DNMT3A mutations were acquired by 8 weeks of gestation to 7.6 years of age in 4 patients, and a PPM1D mutation was acquired by 5.8 years of age. Additional genomic events occurred before or following JAK2V617F acquisition and as independent clonal expansions. Sequential driver mutation acquisition was separated by decades across life, often outcompeting ancestral clones. The mean latency between JAK2V617F acquisition and diagnosis was 30 years (range 11-54 years). Estimated historical rates of clonal expansion varied substantially (3% to 190% per year), increased with additional driver mutations, and predicted latency to diagnosis. Our study suggests that early driver mutation acquisition and life-long growth and evolution underlie adult myeloproliferative neoplasms, raising opportunities for earlier intervention and a new model for cancer development.

Caught in the antiviral crossfire: Ganciclovir-associated mutagenesis in HSC transplant recipients.

In this issue of Cell Stem Cell, de Kanter et al. (2021) show that most allogeneic hematopoietic stem cells do not acquire additional somatic mutations following transplantation. However, they observe somatic mutagenesis associated with the antiviral drug ganciclovir and find plausible evidence that it may contribute to some post-transplant malignancies.

Unmet clinical needs in the management of CALR-mutated essential thrombocythaemia: a consensus-based proposal from the European LeukemiaNet.

Recommendations regarding management of essential thrombocythaemia rely on studies done before the discovery of the CALR mutation. On May 20, 2020, the European LeukemiaNet annual meeting was held with the goal to identify unmet clinical needs in myeloproliferative neoplasms. Because patients with a CALR mutation have specific clinical characteristics, treatment of CALR-mutated essential thrombocythaemia was considered an unmet clinical need by the European LeukemiaNet. The elaboration of a consensus document with recommendations according to current evidence was proposed as a solution for resolving uncertainties in the treatment of CALR-mutated essential thrombocythaemia. A steering committee comprising four European LeukemiaNet members was then formed and a panel of ten experts in the field was recruited. The experts proposed 51 potential unmet clinical needs in the management of CALR-mutated essential thrombocythaemia and were asked to score the relevance of each topic. Those topics that obtained the highest scores as relevant unmet clinical needs were identified, including antiplatelet therapy in patients at low risk, definition of extreme thrombocytosis and its management in patients at low risk, indications of cytoreduction and targets of therapy, first-line treatment of choice in young patients (<60 years), and management of pregnancy. After the steering committee revised the available evidence for each topic, a consensus on management and proposal for improving knowledge was achieved by use of an email-based, two round, Delphi approach. Consensus was achieved when 90% of the panellists agreed with a statement and included 14 recommendations and six solution proposals. Key recommendations included careful observation for asymptomatic patients with classical, low-risk, CALR-mutated essential thrombocythaemia without cardiovascular risk factors; caution in the use of antiplatelet therapy for symptomatic patients at low risk with platelet counts of 1000-1500 × 109 platelets per L, in such cases cytoreduction is an adequate option, especially if adquired Von Willebrand disease is present; cytoreduction is recommended for extreme thrombocytosis (platelet count >1500 × 109 platelets per L) with pegylated interferon alfa being the preferred option for younger patients; both hydroxycarbamide and anagrelide might be given to patients ineligible for pegylated interferon alfa; and treatment algorithms for patients with high-risk pregnancies should not be changed according to genotype. The European LeukemiaNet proposes to use these recommendations in the routine management of patients with CALR-mutated essential thrombocythaemia, and designing new clinical studies in this field might be useful.

Lineage tracing of human development through somatic mutations.

The ontogeny of the human haematopoietic system during fetal development has previously been characterized mainly through careful microscopic observations1. Here we reconstruct a phylogenetic tree of blood development using whole-genome sequencing of 511 single-cell-derived haematopoietic colonies from healthy human fetuses at 8 and 18 weeks after conception, coupled with deep targeted sequencing of tissues of known embryonic origin. We found that, in healthy fetuses, individual haematopoietic progenitors acquire tens of somatic mutations by 18 weeks after conception. We used these mutations as barcodes and timed the divergence of embryonic and extra-embryonic tissues during development, and estimated the number of blood antecedents at different stages of embryonic development. Our data support a hypoblast origin of the extra-embryonic mesoderm and primitive blood in humans.